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Abstract 

The design of iterative schemes for sparse matrix computations 
often leads to constrained polynomial approximation problems on 
sets in the complex plane. For the case of ellipses, we introduce a 
new class of complex polynomials which are in general very good 
approximations to the best polynomials and even optimal in most 
cases. 
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§1. Introduction 

We consider complex Chebyshev approximation problems of the type 


E n(r,c)= min ||p||f r , 

pen„:p(c)=i 


Sr ■= max|p(s)|. 


( 1 ) 


Here n n denotes the space of all complex polynomials of degree at most n, £ r 
is any ellipse with foci ±1 and semi-axes (r ± r -1 )/ 2, r > 1, and c E C \ £ r . It 
will be convenient to express c as a point on the boundary D£r of the ellipse 
£r, R > r, i.e. c = c(R, 7) = ((R + R~ l ) cos(y) + i (R - R~ 1 )sin(j))/ 2 , 7 e 
[ 0 , 2tc). Since Haar’s condition is satisfied, there always exists a unique optimal 
polynomial p„(z;r,c) of (1). 

Problems ( 1 ), in general with £ C C any compact set instead of £ r , arise 
in numerical linear algebra. E.g. the design of iterative methods for the solu- 
tion of large sparse non-Hermitian linear systems Ax = h with best possible 
convergence rates [2], the computation of optimal polynomial preconditioners 
for conjugate gradient type algorithms for Ax = b [ 7 ], or the acceleration of 
eigenvalue methods for A [6] all lead to problems of this type. However, for 
arbitrary sets £ the optimal polynomials are in general not known explicitly 
and therefore the methods are usually based on polynomials which are only 
asymptotically optimal. A popular choice for the set £ are ellipses, and then 
the scaled Chebyshev polynomials t n (z,c) T n (z)/T n (c) are used as approx- 
imations to the optimal polynomials of (1) [ 4 , 6]. Clayton [l] showed that even 
t n (z;c) = p n (z]r,c ) if c is real, and in general t n is nearly optimal for (1) as 
long as n is large. However, in some of the applications we mentioned, polyno- 
mials with small degree are used and typically the distance between c and £ r 
is small. Depending on the position of c on 9£r, ||/„(z; c)||f r > 1 can occur, 
and then t n yields no useful approximation (cf. Example 1 given below). 

In this note, we introduce a new class of asymptotically optimal polyno- 
mials <7„ for Problem (1) which always satisfy ||<7n(~'> c)||*- r < ||t„(z;c)||f r and 
\\q n (z] c)\\ £r < 1. Moreover, they are even optimal in most cases. 

§2. Results 


The q n are defined by 


qn(z]c) 


Tnjz) + On _ 0 . sin(n7) 

Tn(c) + a n ’ an ~-\R”-R-»)- 


(2) 


Here a n is the solution of the extremal problem 


M n (r, c) = min max 
a eC *£S r 


Tnjz) + Of 
T„(c) + o 


( 3 ) 


We summarize the important properties of q n (z,c ) in the following 
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Theorem 1. [3] 

(a) q n {z\ c ) has precisely 2 n extremal points zj, j = 1,2,..., 2n, on <9<? r with 
HtfnO; c)|| fr = M n (r, c) = (r n + r~ n )/(R n + R~ n ). 

(b) There exists a number R 0 (n,r ) such that q n (z\c) = p n (z\r,c) for all 
c 6 with 72 > i?o(^,r). 

('cj Let c € be such that R > r(9r 4 — l)/(r 4 — 1). Then, there exists an 
integer n 0 (r, R) such that q n {z\ c) = />„(■*> r > r ) hor a 11 n > n 0 (r, /?). 

Discussion: Supported by numerical tests (c.f. Example 2), we conjecture 
that (c) is true for arbitrary R. > r > 1. ||9n( z ; f! )||£ r does not depend on the 
position of c on d£ji and ||fl n (z; c)|| fr < ||t n (r;c)|| fr , where equality holds iff 
sin(n 7 ) = 0, e.g. for c 6 IR (cf. Example 1). The proof of Theorem 1 is based 
on the following characterization [5]: q n ( z 'i c ) = p,,(z;r, c) iff the linear system 


2 n 

^^jQn(zj]c)(zj -c)p(zj) = 0 for all P G !!„_! (4) 

i= i 


has a nontxivial and nonnegative solution. See [3] for the explicit solution of 

( 4 ). 

Example 1 . We compare \\q n ( z ] c )\\s r (continous curve) and \\t n (z] c )\\s r 
(dashed curve) where r = 1 . 1 , R = 1.2 for 7 G [ 0 , 7 r] and n = 3, 4 (cf. 
Figure 1 ). 
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Figure 1 . Maximum norm of q n and t u 


The nontrivial solutions of (4) always lead to a lower bound for the min- 
imal deviation of problem ( 1 ), which is sharp in a certain sense: 

Theorem 2. Let 0 j, j = 1,2, ... , 2 n, be any nontrivial real solution of (4), 
normalized such that I a j\ ~ 1 > then 


2 n 


L n (r,c) = 


M n (r,c) 




< E» ('■•<’), 
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where equality holds iff q n {z\c) = p n (~;r,c). 
Proof: Let p 6 II n _j. From (4) we obtain 


2n 


X^i3n(*i;c) = - (Zj -c)p(Zj)) 

> =1 . 7=1 ;=1 

2 n 

< (M„(r,c)^|„ J |) ||l-(r- c )p(z)|| £ ., 


]=t 


and the result follows. ■ 

We illustrate that q n is in general nearly optimal in the following 

Example 2. In this example we compute the relative deviation D n (r,c) — 
( M n (r,c ) - L n (r,c))/M n (r, c) where r = 2 for c E {-2.1, 2.1] x [ — i2.1, i2.l] 
(he re D n (r,c) := 0 if c E £ r ) and n = 2, 3, 4, 5. Note that D n (r,c ) = 0 if 
q n (r,c) is optimal (cf. Figure 2). We obtain max cgC D n ( 2, c) < 0.1024, 0.0498, 
0.0336, 0.0210 for n = 2, 3, 4, 5 resp. 



r = 2, n = 2 
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r = 2, n = 3 



r = 2, n = 4 
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r = 2, n = 5 


Figure 2. Relative deviation of q n 
References 

1. Clayton, A. J., Further results on polynomials having least maximum 
modulus over an ellipse in the complex plane, U.K.A.E.A. Memorandum 
- AEEW - M 348, 1963. 

2. Eiermann, M., W. Niethammer, and R. S. Varga, A study of semiiterative 
methods for nonsymmetric systems of linear equations, Numer. Math. 47 
(1985), 505-533. 

3. Fischer, B. and R. Freund, On the constrained Chebyshev approximation 
Problem on ellipses, J. Approx. Th. , to appear. 

4. ManteufFel, T. A., The Tchebychev iteration for nonsymmetric linear 
systems, Numer. Math. 28 (1977), 307-327. 

5. Rivlin, T. J. and H. S. Shapiro, A unified approach to certain problems of 
approximation and minimization, J. Soc. Indust. Appl. Math. 9 (1961), 
670-699. 

6. Saad, Y., Chebyshev acceleration techniques for solving nonsymmetric 
eigenvalue problems, Math. Comp. 42 (1984), 567-588. 

7. Saad, Y., Krylov subspace methods on supercomputers, SIAM J. Sci. 
Stat. Comp., to appear. 



